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Abstract 

Effect of a "random anisotropy" type disorder on orientation of 
magnetization in an amorphous magnet is considered. It is shown, 
that principal corrections to the free energy of the magnet originate 
from the fluctuation of the order parameter in the directions of its 
degeneracy. The types of disorder are found, which lift the continuous 
degeneracy and suppress the Larkin-Imry-Ma mechanism of disruption 
of the long-range order in continuously degenerate systems. 

1. Continuously degenerate ordered systems react strongly even on a weak 
disorder. It has been demonstrated by Larkin [lj in connection with the 
translationally degenerate system and by Imry and Ma [2J for systems with 
rotational degeneracy that in the spaces with dimensionality d < 4 an arbi- 
trary weak disorder disrupts the long-range order. Let magnetization M be 
the order parameter, £- the correlation length and g-a, small parameter, char- 
acterizing the weakness of the disorder (it will be specified later). According 
to the Imry and Ma argument orientation of M can vary in space to adjust to 
the fluctuations of the disorder on a scale L v ~ £/ ' g. This nonuniform state 
has lower energy then the uniform for a relative value of the order of g 2 . It 
is assumed in this argument that the continuous degeneracy is not lifted by 
corrections to the energy of a lower order in g. The known examples [31 HJ [5] 
demonstrate that some types of disorder can orient M, i.e. to lift a contin- 
uous degeneracy. It was demonstrated also [B] that for the superfluid 3 He in 
aerogel corrections of the lowest order in g can suppress the Larkin-Imry-Ma 
effect. 

In a present paper the argument of Ref . [6] is applied to a more simple ob- 
ject - amorphous ferromagnet with the disorder of the "random anisotropy" 
type [7] in a vicinity of the Curie temperature T c . For such ferromagnet 
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the particular types of disorder are found, which lift the continuous degen- 
eracy and suppress the effect of Larkin-Imry-Ma. As an example the three- 
dimensional ferromagnet is considered with the weak random anisotropy in 
in the x, ?/-plane 0. It is shown that such disorder orients magnetization in z- 
direction. The orienting effect manifests itself in a region where the disorder 
is still small and can be treated as a perturbation [S] . 

2. The Landau expansion of the free energy of an amorphous ferromagnet 
in the powers of M in a vicinity of T c has the form: 

F f = F n + J d"r[TMjMj + b -{MjMjfHl^^+r 1 j l {v)MjM l + Kj l MjM l \. 

(1) 

The coefficients here have the usual meaning: r = (T— T c ) /T c ; b = const., b > 
0, the spatial rigidity is denoted as £q. The units for M are chosen to yield 
the proper dimensionality to the integrand. Random anisotropy is described 
by the term Vji(r)MjMi, where t]ji(r)- a real symmetric random tensor, its 
trace r) nn (r)/3 describes the local shift of the transition temperature and the 
remaining part % (r) = r]ji(r) — 5ji(rj nn (r)/3) - the local splitting of T c for 
different components of M. The ensemble of tensors ^(r) is assumed to be 
spatially uniform on the average. The average anisotropy Kji is written as a 
separate term, i.e. < r]ji(r) >= 0. The results will be expressed in terms of 
the Fourier transforms of binary correlation functions ^(r) in a limit k — > 0: 
$ji mn (0) =< T)ji(k)r] mn (—\c) >|fe=o- For the moment no further assumptions 
are made about these functions. The equilibrium M(r) is found from the 
equation: 

d 2 M- 

tM ] + Kjl Mi + Vji(r)M l - eo-^f + bM 2 Mj = 0. (2) 

To separate the effect of the disorder on orientation of the order parameter 
the = is considered first. Treating rjji(r) as a small perturbation we 

assume that the solution of Eq.(2) has a form: Mj(r) = Mj + rrij(r), where 
rrij(r) is a small fluctuation: \ mj(r) |~| n | 7 M, 7 > 0, its average over the 
volume with a linear dimensions smaller then L v vanishes: < mj(r) >= 0. 
The expansion of Eq. (2) over rjji and mj and separation of smoothly and 
fast varying terms renders equations for m,j(r) and Mj. Possible dependence 
of Mj on r would effect terms of the higher order in rjji then considered here, 
so we can assume Mj = const.. For brevity the bar over Mj in what follows 
is dropped. 

rMj + b[M 2 Mj + 2 < nijmi > Mi+ < mm > Mj}+ < rjjimi >= 0, (3) 
1 the example is suggested by E. I. Kats 
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Tmj + b[M rrij + 2M j M l m l + (mimfimj] - £o^/ = (4) 

The local fluctuations of magnetization m,j(r) in directions parallel and per- 
pendicular to M have to be treated differently. Projection of Eq. (4) on the 
unit vector ft = M/M renders: 

d 2 Tti 

™, + 3&M 2 m„ - il-^r = -V^M, (5) 

where m M = fiimi and rj^ = rjjifijfii. Here and in what follows the summation 
over the repeated Greek indices is not assumed. 

The cubic in mi term can be omitted here. The equation is linear and 
can be solved by Fourier transformation: 

Then for the average fluctuation we have: 

<^»- m, / < ss; > & < 7 > 

The integral here converges at k 2 ~ M/Co ~ V[£CO] 2 - F° r that region we 
assume $ W)1 (k) w $ ww (0), then 



1 gw(0) M 2 
Svr Co 3 ^ 



< m M (r)m M (r) >= _L ^M^ - ? (8) 



For the validity of the above argument the average < m M (r)m At (r) > has to 

be small in comparison with M 2 , or g = _L $ ww(°) <^ i The introduced 

71 ?6v 2 l r l 



here g serves as a small parameter of the present theory. 

To find the transverse with respect to M components of rrij(r) the pro- 
jection of Eq. (4) on the unit vectors A and z>, forming together with ft 
orthogonal basis are used: 

rm x + b[M 2 m x + (m^)m A ] + V^M - tl^r = °> ( 9 ) 

and an analogous equation obtained by substitution of A by v. If the anhar- 
monic term here is dropped the average square of the transverse fluctuation 
would be proportional to the diverging integral: 

< m>( , N ^ tf -/ <y> A, (10) 



The same applies for < m v {r)m u (r) >. To make the expressions for the 
transverse fluctuations finite the anharmonic terms have to be taken into 
account. That can be done in the mean field approximation. In Eq. (9) the 
substitution is made: 

(mimi)m\ ~< m\mi > m\ + 2 < m\mi > mi. (11) 

The sum r + bM 2 also has to be expressed in terms of rrij and rjji with the 
accuracy up to the second order: 

T + bM 2 = -— < rjftimi > -2b < m^m^ > -b < mimi > . (12) 

The average < rj^rrii > is expressed as the integral, diverging at large wave 
vectors k. The diverging part can be included in r and Kji as it was done in 
Ref.pB]. The remaining part of the integral is not singular at small k. Keeping 
only singular terms we arrive at: 

r + b(l + <j x\ + o vv )M 2 = 0, (13) 

where M 2 axx =< rrixmx > and M 2 o vv =< m v m v >. The remaining freedom 
in the orientation of the vectors A and v can be used for turning to zero the 
off-diagonal component Oxv Then Eq. (9) reads as: 

2bM 2 a X xmx - = -M Vx ,. (14) 
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An analogous equation can be written for m u . Eq. (14) has solution: 

Substitution of Eq. (15) in the definition of axx renders the self-consistent 
equation: 

f < T/A M (k)?7v(-k) > d3ft nfi x 
axX J [2bM 2 a X x + e k 2 } 2 (2vr)3' 1 ° j 



which has solution: 



and the same for 



2/3 

. (17) 



2/3 

a w = | ^7^- 1 . (17' 



^U^lUfl (0) 
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The correlation functions $a^a^(0) and $^^(0) have the same order of mag- 
nitude as $ W)J (0). As a result, at g <C 1 in the mean-field approximation 
the transverse fluctuations <taa and of the order g 2//3 , i.e. more impor- 

tant then the longitudinal fluctuation a^. The integral in the self-consistent 
equation (16) converges at ^k 2 ~ bM 2 a\\, or k ~ (l/£(T))g 1//3 , i.e. at g 1 ^ 3 
times smaller wave vectors then those responsible for the longitudinal fluctu- 
ations. The essential wave vectors are still much greater then 1/L V ~ g/£(T). 
This justifies the neglect in the above argument of possible variation of the 
orientation of M resulting from the Larkin-Imry-Ma effect. For g <C 1 the 
conditions of applicability of the perturbative approach cxaa «C 1 and o vv <C 1 
are met. Eq. (13) can be considered as an extremum over M 2 of the effec- 
tive free energy density, which includes the principal corrections due to the 
fluctuations: 

f eff = TM 2 + h(l + (TAA + o vu )M\ (18) 

The inclusion of the corrections results in the substitution of 6(1 + <jaa + ow) 
instead of 6. The sum £ = oaa + o vv may depend on a direction of M. 
Another condition of extremum of f e ff is turning to zero of the derivative 
of £ over the angles, specifying orientation of M. The energy gain at the 
transition is 

26(1 + (TAA + Ovv) 

Corrections <taa and a vv are non-negative. The gain is maximum when both 
oaa = and o vv = 0. That raises a question, for which ensemble rjji(r) 
the maximum gain can be reached. Eqns. (17) and (17') show, that oaa 
and a vv vanish if r}^ = and r\ vll = 0, i.e. if all ^z(r) have at least one 
principal direction in common and magnetization is parallel to this direc- 
tion: f]ji(r)Mi = e(r)Mj. Here e(r) is a random eigenvalue. If all tensors 
r]ji(r) have three principal direction in common the magnetization M can 
be oriented along either of these directions. According to Eqns. (17), (17'), 
(19) the proper orientation of M increases the energy gain for an amount 
~ g 2 / 3 (r 2 /26). This is much greater then the gain at the disruption of 
the long-range order according to the Larkin-Imry-Ma mechanism, which 
is ~ g 2 (r 2 /2b). The lifting of degeneracy suppresses the Larkin-Imry-Ma 
mechanism and the long-range order is preserved. 

The described above lifting of continuous degeneracy by a random pertur- 
bation is analogous to the lifting of degeneracy of the mechanical pendulum 
with the vibrating suspension point (Kapitza pendulum) [S]. If the gravity 
is disregarded the equilibrium position of a pendulum is degenerate. The 
vibrations lift the degeneracy and orient the pendulum so that there is no 
torque in the direction of degeneracy of the pendulum. 
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For the magnet the perturbation has continuous spectrum of the wave 
vectors and it is useful to elucidate the effect of different parts of the spectrum 
on the order parameter. Small and large k in comparison with ~ 1/£(T) 
have to be distinguished. The short wavelength perturbations (k ^> kg) 
effect equally all degrees of freedom. The net effect of these perturbations 
is a shift and, possibly, a splitting of the T c , which can be added to the 
corresponding terms in the free energy. Difference between the longitudinal 
and transverse (Goldstone) fluctuations occurs in the short wavelength region 
(k <C k^). For the transverse degrees of freedom in the three-dimensional case 
the main contribution to the free energy comes from the fluctuations with 
k ~ g 1 / 3 /^. This contribution is g^ 1 ^ 3 times greater then for the longitudinal 
degree of freedom. The disruption of the long-range order according to the 
Larkin-Imry-Ma mechanism also comes via the transverse fluctuations but 
with even smaller k ~ gk^. If the contribution of the transverse fluctuations 
to the free energy depends on the orientation of the order parameter then the 
minimization of this contribution favors the orientation of M, which decreases 
the disrupting effect of the fluctuations on the long-range order. If for some 
orientations of M the singular contribution of the Goldstone fluctuations 
can be eliminated, then the Larkin-Imry-Ma mechanism of disruption of the 
long-range order is fully suppressed. 

3. Consider in more details the case of one common principal direction 
assuming it as z-axis. The eigenvalues for that direction r/ 33 (r) are assumed 
to be statistically independent with the other finite components of ^(r), 
which form a two-dimensional symmetric random tensor i] pq (r), where p and 
q take values 1 and 2. With regard to the i] pq {r) it is assumed that they form 
isotropic ensemble in the (xy)-plane. Their correlation functions at k — > 
can be expressed in terms of two constants: 

% qrs (0) =< rj pq (k)rj rs (-k) >\ k=0 = ® {1) 5 pq 5 rs + & 2) (5 pr 5 qs + ) (20). 

The coefficients $( 2 ) are related to the averages =< ?7ii(k)?722(-k) >|fc=o 
,$ (2) =< J7i2(k)77i 2 (-k) >| fc=0 ,$ (1) +2$( 2 ) =< ?7ii(k)77n(-k) >| fc=0 =< V22 ( k) 7722 (-k) >| fc=0 . 
It follows from these relations, that $ (2) > 0, + 2$ (2) > 0. For sub- 
stitution in Eqns. (18) and (19) the correlation functions have to be ex- 
pressed in the basis A, p,, v. The off-diagonal component <j\ v turns to zero 
if A || fi x z and v — A x fi. Let 6 be the angle between p, and z- axis, 
then < ?7 AM (k)?7 A/1 (-k) >| fe=0 = $ (2) sin 2 #, < ^(k)r^(-k) >| fe=0 = ($ (1) + 
2 $(2) + $0)) sin 2 q cos 2 ? where $(3) =< ^(k)^ M (-k) >| fe=0 $( 3 ) > 0, for a 
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two-dimensional disorder $® = 0. As a result, 

£= ( ) —\ [(^ 2 hm 2 9) 2 / 3 + ((^+2^ + ^) S m 2 9cos 2 9) 2 / 3 } 

(21). 

X = at 9 = as it should be, but 9 = it/ 2 is also extremum. 

This example makes possible to compare the orienting effect of fluctua- 
tions with the average anisotropy, described by the term KjiMjMi in the free 
energy (1). Assume that anisotropy is also uniaxial with the axis, parallel to 
indirection, then Kji is specified by one parameter k: K33 = 2k, k\\ = K22 = 
—k. The trace /% is included in the definition of r. The anisotropy splits T c 
in two - for M in the plane (x, y) and perpendicular to it. The magnitude of 
splitting is At = 3k. Depending on the sign of k one or the other orientation 
of M corresponds to higher T c . Since the random anisotropy always orients 
M along z the case k > 0, when higher T c corresponds to the orientation of 
M in the (x, y)-plane is more interesting. With the account of the anisotropy 

f eff = (r + 2k)M 2 + (r - k)(M 2 + Mf) + h(l + a xx + a vv )M\ (22) 

At (r — k) < f e ff has nontrivial extremum M\ = (M\ + Mf) = — (r — 
K,)/b(l+<j\\), M 2 = 0. It corresponds to the state with the energy f± = — (r— 
k) 2 /2b(l + axx). This state is degenerate with respect to the direction of M in 
the (x, y)-plane and has to have according to Ref. [10] quasi-long-range order. 
At (r + 2ft) < there appears another extremum Mf = — (r + 2k) /b; M\ = 
(Mf + M|) =0. It corresponds to the state with the long-range order and 
with the energy f\\ = — (r + 2k) 2 /2b. This state is more advantageous at 
t pa —6k/<j\\. Transition between the two states according to Eq. (21) is of 
the first order. The requirement of applicability of the used approximations 
imposes rather severe restriction on the anisotropy: 6k/<t\\ 1. 

Consider also the interplay of the disorder with the uniaxial anisotropy 
in the fourth-order terms in M, when without the disorder the free energy 
has a form: 

/ = rM 2 + i[6 || M 3 4 + 6 ± Ml]. (23) 

At r < the energy (23) has two extrema: Mf = —r/bn,M± = and 
Mf = 0, Mf = —r/b±. The more advantageous is the state with the smaller 
b. The disorder effects differently the two states. More interesting is the case 
b± < b\\, when without the disorder orientation of M in the (x, y)-plane is 
favored. The inclusion of the principal order corrections due to the random 
anisotropy amounts to the substitution of 6j_(l + <j\\) instead of b± in the 
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solution. Another solution M || z in the same approximation does not have 
corrections and for &\\ > (fey — it became more advantageous. In the 

temperature interval from r = till the temperature when ctaa = (b\\ — b±)/b± 
the disorder interchanges two states. In the vicinity of T c the state with the 
long-range order is stabilized. Since it has been assumed, that <j\\ <C 1, the 
interchange is possible only for the small anisotropy, i.e. when (b\\ — &_i_)/frj_ <C 
1. 
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